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We present a niechanism to transfer the spectrum of perturbations in a scalar isocurvature 
field ^ onto the matter content in the radiation era via modulated, instant preheating after 
ekpyrosis. In this setup, ^ determines the coupling constant relevant for the decay of a 
preheat matter field into Fermions. The resulting power spectrum is scale invariant if ^ 
remains close to a scaling solution in new ekpyrotic models of the universe; by construction 
the spectrum is independent of the detailed physics near the bounce. The process differs 
from the curvaton mechanism, which has been used recently to revive the ekpyrotic scenario, 
in that no peculiar behavior of ^ shortly before or during the bounce is needed. In addition, 
a concrete and efficient realization of reheating after ekpyrosis is provided; this mechanism 
is not tied to ekpyrotic models, but could equally well be used in other setups, for instance 
inflationary ones. We estimate non-Gaussianities and find no additional contributions in the 
most simple realizations, in contrast to models using the curvaton mechanism. 
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I. INTRODUCTION 

The ekpyrotic scenario is an attempt to provide an alternative to the inflationary 

model of the early universe, see J] for a review. It is based on the moduli space approximation of 
heterotic M-theory in five dimensions. In the proposal, the fifth dimension is sandwiched between 
two orbifold fixed planes and undergoes a bounce, corresponding to a bounce of the scale factor in 
the four dimensional effective description as well. This can provide an alternative solution to the 
horizon problem ^. However, the model is restrained by the singular nature of the bounce, which 
renders the resulting scalar spectral index for the radion ambiguous , even from the five 

dimensional point of view 



Recently, the model has been revived 



11 



12l | by considering additional degrees of freedom, 



such as the volume modulus of the Calabi-Yau. If the fields start out near an (unstable) scaling 
solution, nearly scale invariant fluctuations are generated in an isocurvature mode during the pre- 
bounce phase. These can be converted to adiabatic ones before the bounce by different means; 



for example, the flelds could be moving away from the mu^ 
field ekpyrotic attractor shortly before the bounce 

however 3]). Alternatively, a reflection of fields on sharp boundaries, which can be argued to occur 



ti field scaling solution towards a single 
isl - naturally requiring fine-tuning (see 



naturally in certain models can cause a similar conversion. The idea of employing 



isocurvature perturbations resembles the curvaton mechanism [201 ]. and has already been proposed 
by Notari and Riotto in [21] shortly after problems with the ekpyrotic scenario surfaced. 

However, there remains a subtlety regarding the bounce itself: if the low energy effective theory 
in four dimensions stays valid throughout, energy conditions must necessarily be violated in order 
for a non-singular bounce to occur. The only known mechanism capable of achieving this without 
introducing fatal instabilities is by means of a ghost condensate 1^, [ill , [l^ . whose connection 
to the underlying string theory is unknown. Further, since the conversion of isocurvature to 
adiabatic perturbations needs to occur towards the end of ekpyrosis, both degrees of freedom 
need to evolve quite drastically near the not well understood bounce of the four dimensional scale 
factor Hence one might be concerned that the nature of the bounce has an impact on the 



^ See [13 ] for an updated account on the proposed resolution of the horizon and flatness problem within ekpyrotic 
proposals, attempting to address objections raised in e.g. [j. 

For instance, if the matching of perturbations is done on a constant energy hypersurface, a deep blue spectrum 
results [i,[7a|. 

^ Note that the bounce mentioned in 0] is this sudden change of the trajectory in field space and not the bounce of 
the four dimensional scale factor which is needed to solve the horizon problem. The latter bounce is still assumed 
to be singular in [^. 



3 



spectrum of sea ' 
e.g 



22 



231, 



24 



ar per turbations, an issue studied extensively in four dimensional toy models, see 



26l | and references therein. 
In this article, we propose a mechanism to imprint the spectrum of isocurvature perturbations 
in a field ^ onto the matter content during reheating, that is during the conversion of some preheat 
matter content originating at the brane collision into relativistic Fermions. This provides the 
possibility to generate scale invariant adiabatic perturbations after the not well understood bounce. 
The conversion is achieved by assuming a sensitivity of the coupling constant between the Fermions 
and a preheat matter field on the value of ^. This determination of coupling constants in terms 
of more fundamental fields is a natural feature of string theory. We take the isocurvature field to 
be fixed throughout the bounce and pre-bounce phase in order to remain faithful to the original 
ekpyrotic proposal. In other words, to touch base with the proposals in 



101, 



11 



, IJJ, llJ] , we stay 



close to the multi-field scaling solution. In order to arrive at a scale invariant spectrum for ^, a close 
proximity to this ill-fated unstable solution in the ekpyrotic phase is necessary. As a prerequisite, 
either initial conditions have to be extremely fine-tuned, or an additional pre-ekpyrotic phase needs 
to be added during which ^ is driven towards the scaling solution Ij] . Our proposal is insensitive to 



details of the bounce by construction, and a conversion of isocurvature to adiabatic perturbations 
before the bounce is not needed. 



The employed mechanism of imprinting 



preheating in the inflationary literature 



271, 



perturbation during rehea ting is known as modu 



281, 



29(1 (see also 



30, 



31 



32 



33 



34 



35 



36 



31 



ated 



but it has not yet been applied to ekpyrotic scenarios (even though briefly mentioned in ^]) or 
instant preheating, which we employ. 

In the ekpyrotic proposal, preheating is due to the collision of boundary branes with two salient 
features: flrst, the event is nearly instantaneous, since it is directly related to the collision of branes; 
once they come close to each other, new light states appear and get produced non-perturbatively. 
Second, the force between the boundary branes vanishes near the collision 

To model these features from a four dimensional point of view, we use the scenario of instant 
preheating [39|] where a scalar preheat matter field x> which can be thought of as one of the 
additional light degrees of freedom that appear near the bounce, gets produced once the adiabatic 
field (p is close to zero ^. As (/> evolves towards larger values again, x-particles become heavy and 
decay into relativistic matter, which we take to be a Fermionic field tp. The time of x-decay will 



"* Note that (f> is decreasing in the contracting phase; this ekpyrotic phase ends around 0^0, the field runs quickly 
to ~oo, turns around (corresponding to the bounce), and emerges with <^ > in the post-bounce phase. 
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depend on the value of the isocurvature field if we introduce a dependence of the coupling constant 



271, 



29| 



on ^. This sensitivity leads to modulated perturbations 

The dependency of the coupling constant between preheat-matter and Fermions on the isocur- 
vature field could in principle be computed if the stringy construction, as well as the relation to the 
standard model of particle physics, were well understood. Unfortunately, this is not the case here 
and this contingency is kept arbitrary. As it turns out, this dependence is not overly constrained 
from a model builders' perspective, as long as it is steep enough to yield the dominant contribution 
to the power spectrum. 

To discriminate the above model from 



a 



10 



111 ], we compute non-Gaussianities (NG); we find 



no additional contributions due to the conversion in the simplest case_of an exponential dependence, 
in contrast to set-ups involving the curvaton mechanism (see e.g. 
the recently proposed new ekpyrotic scenarios 



11 



401 1 and references therein) or 



14( 1 . However, larger Non-Gaussianities 



can occur if the dependence of the coupling on ^ is more involved. 

II. MODULATED PERTURBATIONS FROM INSTANT PREHEATING AFTER 

EKPYROSIS 



In single field ekpyrotic scenarios, the radion describes the inter-brane distance from a four 
dimensional point of view. However, other degrees of freedom such as the Calabi-Yau volume 
modulus are present and also expected to be dynamic during the contracting phase (corresponding 
to a decreasing radion). In such a multi- field solution, one can always identify a single degree 
of freedom along the trajectory in field space, the adiabatic mode (j), and fields perpendicular to 
the trajectory, the isocurvature modes. For simplicity, we consider two fields only and denote the 
isocurvature field wit : 



phase as in 



Id, 



11 



1 £, see Appendix I A 21 for details. Given a scaling solution in the pre-bounce 
, [ij, [3], the isocurvature field acquires a nearly scale invariant spectrum of 



fluctuations, as reviewed in AppendixO If the scaling solution is abandoned before the bounce and 
the old, single field ekpyrotic attractor is approached, a conversion of isocurvature perturbations to 
adiabatic ones takes place simultaneously. As a consequence, large non-Gaussianities are produced 
14l ]. To avoid large non-Gaussianities we consider the absence of any conversion before or during 
the bounce. This is only the case if the trajectory in field space is not curved, that is if the 
fields remain close enough to the (unstable) scaling solution throughout Q, [l^. Thus, we take 



^ If the conversion is due to a reflection of fields on sharp boundaries, the resulting NG can be much smaller too 
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^ as approximately constant by construction ^, not because it is stabilized; indeed, a tachyonic 
direction is unavoidable in the new ekpyrotic scenario, if the isocurvature mode is to have a nearly 
scale invariant spectrum of fluctuations [411], see also Appendix lA 21 



Consider now the approach of the boundary branes: we expect the force between them to vanish 
when they are close to each other [3], so that the potential for the radion disappears. Thus, we 
take V{cj},S^) — > for </) ^ 0, in accord with ekpyrotic proposals Q]. After the collision ((/> — oo 
and back again) the universe emerges with an increasing cf). 

Given this generic evolution of ekpyrotic models, we would like to focus on the preheating 
stag e occurring when (/> ~ 0, that is, when branes are within a string length of each other (See 



42 



43 



44 



45l | as well as the review 



46l | and references therein for a sample of the extensive 



literature on (p)reheating). Being in such close proximity, new light degrees of freedom are expected 



to be present 



47|]. For simplicity, to model these new degrees of freedom near (/> ~ 0, we consider 
a single bosonic pre-heat matter field x coupled to the adiabatic field (p. Additionally, to achieve 
reheating of ordinary matter, we couple x to a fermionic field tp. Thus, we take the effective 
potential 

Kff=^9W + /iX#, (1) 

to describe preheating and reheating, for (j) > and increasing. Assuming an appropriate value of 
g, the adiabatic field can produce x-P^-i'ticles when it is close to zero. This production can occur 
rather quickly via non-perturbative effects in a small region around the origin [39], say for </> < cpQ. 
After particle production {(f) > (po) the x particles become non-relativistic. Further, since the 
potential vanishes for ~ 0, the adiabatic field is effectively rolling freely in a universe dominated 
by kinetic energy. Here, we made the conservative assumption that only a small fraction of the 
energy in the adiabatic field gets transferred to x particles initially. Since the equation of state 
parameter for such a stiff fluid is = 1 we have a ~ t^^^. With this time dependence in mind, we 
can integrate the Klein-Gordon equation for (p to 

cp=-^ln{t/to) + (Po, (2) 
zv ovr 



where to = I/SHq ~ Mp/2v37r(^o is the Hubble time around particle creation, depending on the 
initial velocity (pQ > oi the adiabatic field just before the production of x-particles. The energy 



^ We make this assumption primarily to focus on the effects due to modulated instant preheating, not because it is 
needed for the proposed mechanism to work; one could conceivably abandon the scaling solution before the bounce 
and still imprint perturbations during preheating as discussed in this article. 
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in the adiabatic field redshifts very fast oc thereafter. Since (po <C Mp, we are justified to work 
in the 'instantaneous' limit, that is with 0o ~ 0; this corresponds to the fast creation of particles 
once the branes are within a string length of each other. We can also neglect backreaction of 
X particles, which would become important at ti when (p ~ Mp, since x particles decay earlier 
for appropriate values of h at tc < ti. Up until tc the energy in the preheat matter field is 
increasing, due to the increase of (p to which it is coupled. If no decay channel were open to x the 
adiabatic field would bounce back due to backreaction, providing an example of moduli trapping 



47|, 



48 



49 



3, 51, 52, 

This is not the case here, since x-pa-rticles can decay into Fermions. 



Thus, the time tc marks the transition to a radiation dominated phase. In the cyclic scenario, 
the interbrane potential for the radion would kick in at some point later, slowing it down until it 
eventually turns around in the far future to commence another cycle j^, [5^. However, it is not 
clear to us how a cyclic evolution occurs in the presence of more fields, especially in the presence 
of an unstable direction. Accordingly, we will be content with a single bounce m this study. 

The preheating stage described above is the instant preheating scenario of [39^ proposed within 
the inflationary framework; we refer the interested reader to [39] for more details, e.g. regarding 
constraints on coupling constants g and h such that sufficient production of x and ip particles takes 
place. The new feature above is the implementation of instant preheating within the ekpyrotic 
framework. 

Next, consider a dependence of the coupling "constant" h on the isocurvature field ^, which has 
been a spectator during the bounce and the preceding ekpyrotic phase ^. Such a dependence is 
indeed expected in many stringy models, but a specific construction is lacking so far. Taking this 
uncertainty as an opportunity for optimism, we regard /i(^) as a free function, constraint by the 
requirement that h{S,scaiing) yields sufficient production of Fermions. That way we can provide a 



271, 



28 



ft 



the decay rate of x-particles 



new realization of modulated perturbations from reheating 

depends on this coupling and the radion field value via F oc h'^(p [39]. Since (p is only increasing 
logarithmically, the decay rate before ti is very well described by its value at ti. Henceforth, the 
only dependence on the fermionic coupling is 

r(xh^, (3) 

so that the time at which Fermions are produced becomes 

UO^HO'^- (4) 



^ In principle g could also depend on ^, but since g is heavily constraint from requiring that x-particles get produced 
sufficiently, we will ignore such a dependence. 
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The corresponding adiabatic field value at tc depends only logarithmically on h via 



2V37r 



ln(/i 



-2^ 



+ ... , 



(5) 



where the dots denote terms independent of h. The dependence on h{^) in ([3D induces an additional 
contribution to the power spectrum and higher order correlation functions. 

To compute these effects, we use the dN-iorma : 
extended by Sasaki and Stewart [s^] and others 



58 



ism, first proposed 



591, 



60, 



61 



62, 



631, 



by 



64 



Starobinsky in 



65!, 



5g] and 



66( 1 ■ This formalism 



relies on the fact that on super-Hubble scales perturbations evolve classically (no oscillations) and 
fluctuations can be taken smooth on scales smaller than the Hubble radius. Based on this, one can 
relate the perturbation of the volume expansion rate 6N to the curvature perturbation (which 



67|, 



Given this 



is conserved on large scales in simple models) even beyond linear order 
relationship, correlation functions of (k denotes a Fourier mode) can be computed in terms 
of the change of during the evolution of the Universe. The formalism cannot only be applied 
to simple models of inflation, but also to scenarios where fluctuations are imprinted on after 
reheating, e.g. the curvaton scenario [6d|, and provides an extremely efficient way of estimating 



non-Gaussianities (here the non linear (5A^-formalism developed in 



59|, 



601 1 is needed). 



Following this thread of thought, the resulting power-spectrum for two fields becomes 



dN 



where the derivatives are to be taken deep inside the pre-bounce phase. Extending the 5N- 
formalism to more fields and higher orders results in 

C{t2,x) = V —5<Pi + - V ——6cl)i64>j + ... , 7 
^ d(j)i 2 ^ d(j)i d(t)j 

so that higher order correlation functions, such as the bi- or trispectrum, can be computed too (see 



dN 



Pa 



(6) 



e.g. 



63 



Consider now the number of e-folds from the production of x-particles until after their decay 
into Fermions, when we are deep in the radiation era, specifically at tr 

Nr{h) = f Hdt 



to 



t-(h) 



H_dt + 



to 



Vc (h) 



drj + 









dt 


it+{h) 




r !± 


drj 




Jr,t{h) V 





(8) 



Note that the separate Universe formalism developed by Rigopoulos and Shellard in e.g. [631 equivalent to the 
(5iV-formalism. 
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Here, we switched to conformal time and used the Hubble parameter TC = a' /a = r/r] with 
r_ = 1/2 (kinetic driven phase) or r-|_ = 1 (radiation phase). Conformal time makes a jump at the 
transition, r]j_ = r/~r_)_/r_ = 2r/~, since Ti. is continuous at the transition via the Israel junction 
conditions Q] • Integrating ([8]) yields 

A^^ = r_ ln{7]-/r]o) - r+ ln(??^" r+Zr/^r.) . (9) 

Taking the derivative with respect to ^ and using r]~ oc tc ^^^^^"^ oc we arrive at 

dNr _ , .1 dri~ dh 

= (r_-r+) — ^ — 

^ , ^ 1 dh , . 

= 2(r_-r+)(r_ + l)-— , (10) 

so that the additional contribution to the spectrum becomes 



d^ 

2(r_-r+)(r_+l)i|^^) 
3 1dh^^ 



It should be noted that the more pedestrian method of using the Deruelle-Mukhanov matching 
conditions across a sharp transition of the equations of state parameter could be used 

alternatively |71|, but the 5A^-formalism is more efficient. The contribution in (jlip can dominate 
over the one of the adiabatic field, if there is a steep dependence of h on the primordial value of 
the isocurvature field Therefore, the relevant spectrum for the cosmic-microwave background 
radiation does not need to be the one of the adiabatic field, but can be provided by an isocurvature 
field (see Appendix I A 21 for the computation of its spectral index) ^. 

The difference to the curvaton scenario [2^ or the recently proposed new-ekpyrotic scenario 



Id, 



111 ) is simple; namely, to employ the curvaton mechanism one needs the additional degree 
of freedom to be dynamic in a very peculiar way close to the end of ekpyrosis /inflation, with the 
danger of being sensitive to the bounce in ekpyrotic models. In our scenario, ^ remains constant 
throughout ekpyrosis, causing the decay of the pre-heat field x fo occur at slightly different times, 
in accordance with ,^'s (isocurvature) fluctuations. Thus, the proposed mechanism is rather robust 
and insensitive to the nature of the bounce by construction. 



During completion of this work, [72l | appeared providing an application of the (5A''-formalism to modulated pertur- 
bations in ordinary preheating after inflation - see note added after the conclusions. 
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The main drawback of the ekpyrotic mechanism to refurbish the isocurvature field ^ with a 
nearly scale invariant spectrum of fluctuations in the contracting phase is the necessity of ^ to have 
a tachyonic direction [41]. As a result, ^ is not stabilized, but it is supposed to linger near the 
crest of the potential. Even if our interest lies in a single bounce only, to be close to the scaling 



solution 



13 



14 



4l| there is a need for highly fine-tuned initial conditions. As long as no compelling 
mechanism is found to explain such initial values, the model remains rather unpleasant, proving, 
once again, the difficulty of constructing alternative mechanisms to inflation with regards to the 
observed spectrum of fluctuations. One can by all means add a pre-ekpyrotic phase and postulate 
a potential such that ^ is driven very close to ^scaling 



121, but it remains to be seen if such a 



potential arises in a concrete stringy construction. 

We would like to emphasize that modulated perturbations from instant preheating as described 
above cannot only occur after an ekpyrotic phase, but they can also be present in an inflationary 
framework; if an isocurvature mode has a scale invariant spectrum of fluctuations after inflation 
(which is the case if additional light degrees of freedom are present during inflation), the desired 
scale invariance of ^ after reheating can be achieved without relying on the curvaton mechanism. 



This is only a slight extension of 



27 



28( 1 Meanwhile one may naturally wonder how the above 



mechanism, which converts isocurvature perturbations to adiabatic ones, could be differentiated 
from a conversion before the bounce or from simple models of inflation. A possibility is to con- 
sider non-Gaussianities, namely, the non-linearity parameter fi\iL characterizing the ratio of the 
bispectrum divided by the square of the power spectrum. The observational bound on the full non- 
linearity parameter obtained from the WMAP3 data set alone is —54 < /jy^ < 114 
future experiments will improve upon it considerably [77 



, I75I, 



74|,l75|,l76|, and 



contribution to fi\fL can be computed within the 5A^-formalism to [64] 



, Q, 0i • The momentum independent 

n 



TJNL 



J2jjNiNjNij 



where the sum runs over all fields, (j) and ^ in our case, and = dN/d^, 
from modulated perturbations after inflation have been considered in the past 
the (5A^-formalism. A recent study 



(12) 

Non-Gaussianities 



80] , without invoking 



8l[] confirmed, in the same framework, the equivalence of results 
based on the (5A^-formalism. Hence, without further ado, we apply the (^A^-formalism to the case 
of instant modulated preheating. 



See also [73(] for a related model, where the decay rate depends on the phase of a complex scalar inflaton field. 
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According to p2|) with pO|) . the additional contribution to /nl due to ^ is proportional to 

^«C = ?(7?-72) , (13) 



where we defined 



1 (9"/i 



7n = r- 



(14) 



If we take a simple exponential dependence h ~ e^^, this additional term vanishes identically. 
Further, it is easily checked that the same holds true for higher order correlation functions, such 
as those measured by the non linearity parameters t^l and gNL of [63|] (see also (81!]). Thus, the 



most simple realization of the above mechanism is not expected to generate non-Gaussianities in 
addition to the intrinsic ones of the isocurvature field, which can still be quite large in ekpyrotic 
models {lol.This has to be contrasted with the negligble non-Gaussianities in single field inflationary 
models (see however [s^ for more intricate single field models), the moderately large NG in the 
curvaton scenario 



Q o. afte. a .efleCion of «elds o„ a sl,a.p boundary i. ekpyrojic .od^s Q, 

and the large NG due to a conversion during the ekpyrotic phase as examined in [11|, [Ij, [15|] . 

As an example of the latter, consider a transition from the unstable two field scaling solution 
(reviewed in Appendix E]) to the single field, old ekpyrotic attractor solution, as discussed in jisl : 
the curved field trajectory responsible for the conversion causes an evolution of perturbations even 
after horizon crossing. This evolution is the primary cause for the generation of non-Gaussianities. 
Again, using the 5A^- formalism, one can compute /nl analytically to [if 



/i = ^cL (15) 

where q is the exponent in the potential ()Aip of the field which becomes subdominant at late 
times Since q has to be reasonably large in the ekpyrotic phase in order for fast roll to occur, 
the non-linearity parameter becomes quite large too; indeed, current observations put the model 
already under pressure [l^. Accordingly, such a conversion before the bounce will be refuted if 
large non-Gaussianities are not found in the near future. 

A conversion due to a reflection of fields on a sharp boundary, during a phase where scalar 
kinetic energy dominates, leads to Jnl ~ C(ci), smaller than (jlSh but still potentially measurable 
in upcoming experiments |l3] . Our proposal with h ~ e^^ falls into the same class. Of course, more 
complicated /i(^) could lead to potentially large non-Gaussianities even in our proposal. Assuming 



This result was also confirmed numerically in [Tst - 
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^ Nff,, we can estimate the additional contribution to the non-Unearity parameter to 

/Si ^ I) . 

where we used and (fTUj) in (fl^ . At this point a better understanding of /i(^) is needed in order 
to compute 7^ and provide more concrete predictions. The same holds true for the amplitude of the 
power spectrum [831]. In absence of such a study, which requires a much improved understanding of 
the underlying string theoretical construction, one might favor the simple exponential dependence 
put forward here with ~ 0. 

Before we conclude, we would like to mention the possibility to discern modulated preheating 
from standard single field inflation by testing the consistency relation between the amplitude of 
curvature perturbations, the amplitude of tensor modes and the tensor spectral index; in modulated 
reheating, this relation need not be the same as in simple inflationary models [32]. Again, a better 
understanding of the underlying stringy construction of the new ekpyrotic scenario is needed in 
order to make any concrete predictions. 



III. CONCLUSION 



In brief, we presented a mechanism for imprinting scalar perturbations onto the matter content 
after the bounce in ekpyrotic models of the universe (or after an inflationary phase), providing 
a realization of modulated perturbations from instant preheating. The spectrum is given by the 
one of a scalar isocurvature field, which determines the magnitude of the coupling constant be- 
tween preheat matter and Fermions. The model is independent of the details of the bounce in 
ekpyrotic models, as well as the spectrum in the adiabatic mode before reheating. We estimate 
non-Gaussianities and find no extra contribution in the most simple realization of modulated in- 
stant preheating, in contrast to models relying on the curvaton mechanism in one way or another. 
This makes our proposal appealing, especially as long as primordial non-Gaussianities are not 
observed. 



Note added: during the completion of this work a related article by Suyama and Yamaguchi appeared in (8l| |, 
discussing non-Gaussianities in models of modulated reheating after inflation with multiple isocurvature fields. 
Their study is based on the ^A'^-formalism, similar to the approach followed in this paper. The difference lies in 
the focus on inflation and the absence of instant preheating. Therefore, the imprint of fluctuation onto the matter 
content in [sjj occurs over many oscillations of the inflaton fleld during ordinary preheating, so that the transition 
to a radiation dominated epoch is less sharp. In fact, this is the origin of the nontrivial Q{T/Ho) in Isil] . which 
complicates the computations and yields potentially large NG. We agree with the approach taken in [SlJ and it 
should be straightforward to generalize their results to ordinary preheating after ekpyrosis. 
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APPENDIX A: GENERATION OF SCALE INVARIANT ISOCURVATURE 
PERTURBATIONS IN A CONTRACTING UNIVERSE 

For pedagogical reasons, we would like to review the origin of scale invariant perturbations in 
an isocurvature mode during a contracting phase of the universe in the most simple cases, following 



closely 



1. Background Dynamics 

Consider two scalar fields i = 1,2, with canonical kinetic terms and uncoupled, exponential 
potentials 

2 

W = -Y^ Vie-^'^P' , (Al) 

i=l 



with Cj = ^j2/pi > and Vi > 0. Their equations of motion are 



ipi + 3Hipi = (A2) 



and the Friedman equations read 



m' = w + j2l^l (A3) 

i=l 

H = -j2l^l (A4) 
i=i ^ 

Here and in the following we set SvrG = 1. Given the potential in (|Aip . there are two attractor 
solutions with either Lp\ or ^p^ dominating. Each attractor corresponds to a single field ekpyrotic 
solution. However, there is an additional, exact, multi-field scaling solution to the equations of 
motion, for which the ratio of the fields' kinetic and potential energy remains constant. This 
solution is easily obtained after a field redefinition into an adiabatic field 

<t> ^ 2^^±^ (A5) 
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and an isocurvature field 



so that the potential becomes 



where = Ci + C2 and 



(A6) 

W = -C/(Oe""^ , (A7) 



U{0 = Vie-^"'/"^^"^ + Fse^'^^/^^^'^^ . (A8) 
One can check that —U has a maximum at 



eo = ^.^ln(-i-^) . (A9) 



The scaling solution corresponds to ^ = .^O) whereas (p rolls down the exponential potential 
Obviously, this solution is unstable, and it is indeed the tachyonic instability which gives rise to 
a scale invariant spectrum of perturbations in ^, given that the universe is contracting with an 
equation of state parameter satisfying w ^ 1. The presence of an instability is not a peculiarity of 
this particular potential, but a necessity [41^]. Naturally, such multi-field ekpyrotic proposals are 
rather contrived, as long as no compelling reason is found to start with initial conditions in close 
proximity to the repeller; one could postulate an additional pre-ekpyrotic phase [l^, to drive ^ 
towards the repeller at early times via an additional potential with positive mass squared and a 
coupling to Fermions Ignoring this problem, let us give the model the benefit of the doubt and 
assume that ^ stays close to ^0 throughout, so that ^ ~ and U{(,) ~ f^(^o) = const. It follows 
that the Friedmann equation yields power law contraction a oc {—ty with p = It will be 

useful to introduce 

e-=^(l + w;) = ^ = ^»l, (AlO) 
where we used the fast roll parameter 

in the last step (here (•)^(^ denotes a partial derivative with respect to (f)). Further, the Klein Gordon 
equation for the adiabatic field yields (j) = Afterwards, we will switch to conformal time 

rj = J 1/ a dt so that a oc {—r]y^^^~^\ 



In case the fields move away from the scahng solution, one can still define instantaneous isocurvature and adiabatic 
fields. Since we assume the validity of the scaling solution up until the end of preheating, the initial decomposition 



remains valid and we refer to (f> and ^ as the adiabatic and isocurvature field throughout. 
The coupling to Fermions is needed in 
would amplify oscillations otherwise. 



^'^ The coupling to Fermions is needed in order to dissipate energy while ^ is oscillating around ^0, because blueshifting 
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2. Perturbations 

A useful quantity is the Sasaki-Mukhanov variable, for multiple fields defined via 

Qi = 5h + '^i>- (A12) 

Here, ip is the Newtonian potential (the diagonal metric perturbation in the spatial directions). Qi 
has the easy interpretation of the I'th field perturbation in the spatially flat gauge (■0 = 0). The 
comoving curvature perturbation IZ is related to the Qi via TZ = 4>iQilYl] Fourier 



components of Qj satisfy (see e.g. 8j] for a derivation from first principles) 



Qi + 3HQj + ^Qj (A13) 



a2 



J 

It is already obvious at this stage that one should generally expect a coupling between the different 
Q/'s. However, working in the = gauge and decomposing perturbations into an adiabatic one 
along the field trajectory and an isocurvature one perpendicular to it, it becomes evident that the 
isocurvature mode seeds the adiabatic mode only if the trajectory in field space is curved 8J|. Since 



our scaling solution corresponds to a straight line in field space, no conversion of the isocurvature 
mode to the adiabatic one will occur and the equations decouple if they are written in terms of 
= 5(f> and = 5S,. 

To be precise, using ^ = as well as W^^ = we arrive at 

= 6^ + 3H6C + + > (A14) 

= 6<P + 3H5^ + + W^^^ - ^ {jj(t>''^ ) ^<l^ ■ (A15) 

Introducing = a6^, u^f^ = a6(/) and working in conformal time, the above equations reduce after 
some algebra to 

= u'l+ (^k'^ - ui (A16) 

with J = X, (/> and 

(.--1)^ ' (^1^) 
' - ' (A18) 



{e-lf 



On large scales, TZ is identical to —C,. 
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where we used the background solution of the previous section. 

Imposing vacuum initial conditions well inside the horizon {k'^rj^ ^ /?/ so that /?/ can be 
neglected) we have 

ui = -^ (h{\L)e-'^'^ + 6+(k)e^^'?) , (A19) 



/2k 

where h satisfies (6(k)6^(— k)) = (5'^(k — k) and 6|0) = 0. If we keep the j3i term in the equation of 



(C(k)i/(1)(-A:r?) + C+ {-\^)HfX-kr^)) , (A20) 



motion we arrive at 



ui 

(1 2) 

where H{,/ are Hankel functions of the first and second kind respectively, with index 

Ui = ^yTTWl. (A21) 

Expanding the Hankel functions for large arguments results in 

^ /Ae±^(-(i+2.,)./2) _ .^22) 

^ \ TTX 

Here, the "+" corresponds to hI]^ and the "— " to Using this expansion, we can match uj 

from (|A20]l to (fAT9]l and determine C to 

C(k) = ./^ef(2.,+i)j(l^)_ (A23) 
V 4k 

After horizon crossing, that is in the limit —krj but — In^—krj) not too large, we can use the 
small argument limit of the Hankel functions 

H(]''\x) ^ ±T{uj) - , (A24) 

^ \X J ITT 

where F is the Gamma function. Bringing everything together, we arrive at 

UI = -^2'''T{ui){-krj)^-'''bik) (A25) 
V47r/c 

where 

6(k) ^ 1 (ex(2-i+i)^(k) - e-x(2-/+i)6+(-k)) (A26) 

satisfies also (6(k)6+(k)) = (^^(k — k). Thus we have 6(j) = u^/a oc k^'^'i' and 5S, = u^^/a oc k^'^^. 
Since the power spectrum is defined via 

97r2 1 

6{k - ^')^ri ^ ^{uj{k)uj{k')) (A27) 
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we can easily read off the scalar spectral index to 

= 3-2ui. (A29) 

If we further expand lyj from ()A2ip with ()A17p and (jAlSP in terms of the fast roll parameter 
e = l/2e <C 1 we arrive at 

~ 2 + 4e , (A30) 

~ 4e . (A31) 

As expected, the adiabatic mode carries the deep blue spectrum from single field ekpyrosis, 
whereas the isocurvature mode has a nearly scale invariant spectrum with a small blue tilt. How- 



ever, if slightly different potentials are considered, the index can easily acquire a red tilt 
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At this point we would like to emphasize again that the tachyonic instability of the scaling 
solution is crucial for obtaining a nearly scale invariant spectrum in a contracting universe; to see 
this, let us re-examine how n^ — 1 ~ arose above: from n^ — 1 = 3 — 2^^ we see that ~ 3/2 has to 
hold, which in turn requires (3^ to be close to 2, via the definition of (jA2T]) . If we now look back 
how arose in the equation of motion ()A16p . we see that the leading order contribution during 
fast roll {Ti and Tl' are small) stems from /?g ~ —W^^^a'^r'^ . As a result, the second derivative 
of W has to be negative, which in turn means that ^ exhibits a tachyonic instability (recall that 
W^^ = 0). This is true more generally |4l|, and not an artifact of the simple model employed here. 
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